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We perform a systematic investigation on an asymmetric zig-zag spin ladder with inter-leg ex- 
change Ji and different exchange integrals J2 ± 5 on both legs. In the weak frustration limit, the 
spin model can be mapped to a revised double frequency sine-Gorden model by using bosonization. 
Renormalization group analysis shows that the Heisenberg critical point flows to an intermediate- 
coupling fixed point with gapless excitations and a vanishing spin velocity. When the frustration is 
large, a spin gap opens and a dimer ground state is realized. Fixing J2 = Ji/2, we find, as a function 
of 5, a continuous manifold of Hamiltonians with dimer product ground states, interpolating between 
the Majumdar- Ghosh and sawtooth spin-chain model. While the ground state is independent of the 
alternating next-nearest-neighbor exchange S, the gap size of excitations is found to decrease with 
increasing S. We also extend our study to a two-dimensional double layer model with an exactly 
known ground state. 



I. INTRODUCTION 

Quantum spin ladder systems have attracted much at- 
tention in the past few years^'^. Strong quantum fluc- 
tuations prevent any long range antiferromagnetic (AF) 
order in quasi-one dimension. The magnetic phases of the 
ladder systems are rich and strongly dependent on their 
geometric structures. Several types of disordered "quan- 
tum spin liquid" phases are known^'^'^. Typical examples 
of two-leg ladders are the railroad ladder and the zigzag 
ladder. The railroad ladder has a singlet ground state 
with elementary triplet excitations (magnons)^. Depend- 
ing on the ratio of the leg to rung exchange integrals, the 
zig-zag ladders may have gapless ground states with al- 
gebraically decaying spin correlations or spontaneously 
broken dimerized ground states^. The gapped dimer 
ground state is degenerate, and the elementary excita- 
tions are pairs of spinous. The spin ladders have been 
studied experimentally in compounds such as SrCu203 
and CuGeOs^. 

The two-leg zigzag ladder, which has been well investi- 
gated, is perhaps the simplest example of the frustrated 
spin model and highlights the role played by frustration. 
However, less attention has been paid on asymmetric spin 
ladders where the exchange integrals on both legs are dif- 
ferent. Only the extreme case where one leg of a zig-zag 
ladder is missing entirely (sawtooth or A-chain) has been 
solved^'*. In this paper, we perform a systematic study 
of an asymmetric zig-zag spin ladder, which is a Heisen- 
berg model defined on the structure shown in Fig. 1. In 
general, it is convenient to represent the zigzag ladder as 
a spin chain with nearest-neighbor (NN) exchange Ji and 
next-nearest-neighbor (NNN) exchange J2 corresponding 
to the inter- leg exchange Ji and intra- leg exchange J2. 
Equivalently, the asymmetric ladder model can be repre- 



sented as a chain with an alternating NNN exchange 

H = Y, {Ji^i ■ Si+i + [J2 + i-iys] Si ■ S1+2} , (1) 
I 

where Ji = 1 and J2 ± (5 are the nearest-neighbor (NN) 
and alternating NNN coupling constants, respectively. 
The introduction of S makes the exchanges on top and 
bottom legs different. (5 = is the ordinary zig-zag 
ladder or frustrated spin chain^ and S = J2 is the ex- 
treme case with one leg completely missing. Fixing 
J2 = ^1/2 gives an exactly solved continuous manifold 
of Hamihonians^. The Majumdar- Ghosh (M-G)i° and 
the sawtooth chain^'* are extreme cases with 6 = and 
(5 = J2 = Ji/2 respectively. 
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FIG. 1. The asymmetric zig-zag spin ladder with different 
exchanges on the top and bottom legs. 

For (5 = the model (1) is well understood for gen- 
eral J2^. Frustration due to J2 is irrelevant when J2 < 
J2c^'^^, and the system renormalizes to the Heisenberg 
fixed point^^, whose ground state is described as a spin 
fiuid or Luttinger liquid with massless spinon excita- 
tions. As J2 > J2C, the frustration term is relevant and 
the ground state is doubly degenerate. Particularly, the 
ground state has simple product form at J2 = 0.5Ji^°. 
The critical value of J2C — 0.2412Ji can be determined 
numerically^'-'^^. For J2 > 0.5Ji, quite different field the- 
ory treatments are required depending on the ratio of 
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Since relatively less is known about the asymmetric 
ladder, we study physical effects brought about by the 
leg-asymmetry. Before solving the quantum problem, we 
start with the classical problem which can give us an in- 
structive insight into properties of the asymmetric spin 
model. In the classical limit, the ground state of this 
model is a Neel state for J2 < Ji/4 and a spiral with a 
pitch a = arccos(— Ji/4 J2) for J2 > Ji/4. Both groimd 
states and the critical ratio of exchange integrals sepa- 
rating them, are independent of S. The excitations may 
depend on S, however. 

Certainly, the quantum case is much more compli- 
cated. When the asymmetric exchange interaction is in- 
troduced, some unexpected phenomena will appear. We 
found that the Heisenberg fixed point is no longer stable 
and flows to an intermediate-coupling fixed point with 
gapless excitations and a vanishing spin velocity.^ How- 
ever, there is still reminiscence of classical results. It was 
found that the ground state is independent of S when 
J2 = -^1/2, but the excited gap is decreased by 6. Part of 
the work on weak frustration regime has been reported in 
our previous letter^ , and here we quantitatively study the 
crossover of the excitation spectrum from the symmetric 
M-G model to the extremely asymmetric sawtooth model 
in detail. An extension to a two-dimensional double layer 
model is also presented. 

The outline of this paper is as follows. In section II, 
the effective low energy theory of the asymmetric ladder 
is derived by using bosonization. We qualitatively dis- 
cuss the effect of alternating NNN interaction and com- 
pare our model with the well-known spin-Peierls model. 
The phase diagram of the system is discussed with the 
help of the renormalization group analysis. In section 
III, we study the asymmetric model at the special point 
J2 = Ji/2. The ground state and excitation properties 
are also discussed. In section IV, we generalize our model 
to a two-dimensional double layer model, whose ground 
state is a dimer product state and the excitations are 
magnons. Section V contains our conclusions. 



II. EFFECTIVE LOW-ENERGY THEORY AND 
RENORMALIZATION GROUP ANALYSIS 



where is a bosonic phase field and 11(2:) its canon- 
ically conjugate momentum. Here a is a short-distance 
cutoff, (73 (X 1 — J2/J2C is the Umklapp-scattering am- 
plitude and 51 oc 5 is the amplitude of the alternating 
NNN field. The parameters u and K are the effective 
spin velocity and coupling constants which are given by 



u 



K = 



27r + ff4 - 92 
27r -h 54 + 52 ' 



In genereal, these values are only valid near the free 
fermion point (K = 1), whereas K = 1/2 is fixed by 
the symmetry at the isotropic point. The correspond- 
ing spin-correlation functions can be calculated from the 
boson representation, which gives 



(4) 
(5) 



It is clear that the SU(2) symmetry is restored at the 
isotropic point with K = 1/2. 

For small 93 and gi. Hi could be considered as a per- 
turbation to Hq. Without the gi term, H = Hq + H\ 
represents a standard sine-Gorden model^. The 53 term 
is either marginally irrelevant, which leads to the weak- 
coupling Heisenberg fixed point, or relevant, which drives 
the system to a strong-coupling dimer state. Therefore, 
what we are interested in is how the alternating NNN 
interaction gi changes the physical properties of the sys- 
tem. Qualitative results on the influence of the new in- 
teraction [gx) can be obtained from scaling analysis and 
physical considerations. 

An important information about perturbative opera- 
tors is whether they are relevant, marginal or irrelevant. 
In general, only the most relevant perturbation is impor- 
tant, because the irrelevant operator will scale to zero 
at large lengths. We can give an approximate estimate 
by comparing scaling dimensions of the given operators. 
It follows that e*"**^^^ has scaling dimension of and 

9x^e*"*(^) has scaling dimension of -t- 1. Therefore, 
the scaling dimensions of the Umklapp and the alternat- 
ing NNN terms 53 and gi are 



Following the general procedure of transforming a spin 
model to an effective model of continuum fleld^^'i^, we 
convert the spin Hamiltonian to a Hamiltonian of spin- 
less fcrmions using Jordan- Wigner transformation, then 
map it to a modified Luttinger model with Umklapp and 
backscattering-type interactions. Using the standard dic- 
tionary of bosonizationi^'i'', we obtain the effective boson 
Hamiltonian H = Hq + H\ with 



Ho= dx 
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do 



4K, 



K+1. 



(6) 



At the isotropic Heisenberg fixed point, gs is marginal 
with rfgg = 2, while the gi-term with dg^ = 3/2 is rel- 
evant. We conclude that gi destabilizes the isotropic 
Heisenberg fixed point and the spin liquid ground state. 

On the other hand, for J2 > J2C, the 513 term is 
marginally relevant and renormalizes to a strong coupling 
fixed point in the long- wavelength limit. Near the strong 
coupling dimer fixed point, the gs term is much more rel- 
evant than the gi term. Usually (e.g. 93 ±00), the 
boson field $(a;) locks into a constant value with small 
fluctuations, and an associated excitation gap. Here, the 
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constant solution <I> = ±7r/4 corresponds to the degen- 
erate ground state at the strong couphng fixed point for 
53 — > oo. The standard cos 4$ sine-Gordon equation has 
a pair of solutions of kink and antikink, which describe 
the elementary excitations (a pair of spinon) for the de- 
generate dimer phase. Even including the less relevant 
gi term, the soliton solutions will survive. However, the 
phase locking of gi term is forbidden by the dx^-pre- 
factor to the cos(2$)-tcrm in Hi. In this sense, there is 
no standard strong coupling theory for the gi-term. 

From known results on the sawtooth chain^'^ and the 
M-G model"'^", we expect that the gi term, induced by the 
alternating NNN interaction, does not confine spinous 
and plays a quite different role than the dimerization 
by other degrees of freedom. Moreover, the difference 
in the size of excitation gaps in these two models implies 
that the gi term quite generally competes with the Umk- 
lapp term whereas an external NN dimerization would 
cooperate. As we will show, it turns out that gi opens 
no spin gap despite being a relevant perturbation of the 
Heisenberg fixed point. This result is also corroborated 
by the absence of a magnetization plateau in our model in 
small magnetic fields"'^*. For an alternating NN exchange, 
a magnetization plateau is observed in small magnetic 
fields, but for alternating NNN exchange, it is only ob- 
served in high fields^^'^^. 

To get an instructive insight, we would like to compare 
the alternating NNN interaction in our model with the 
alternating NN interaction in the well-known spin-Peierls 
(SP) model. In the language of field theory, the exter- 
nal dimerization corresponds to a relevant term Hgp oc 
J dx{gsp sin 2$). This term is always much more relevant 
than the 53 term and lifts the degeneracy ($ = ±7r/4) of 
the ground state. With the presence of the g^p term, the 
lowest energy configuration is — 7r/4 if gsp > or 7r/4 if 
gsp < 0. Since the most relevant term is sin 2$, the cor- 
responding sine-Gordon system has a pair of soliton and 
antisoliton solutions (excitations with Sz = ±1) as well 
as two breather solutions (excitations with Sz = 0).^"'^^ 
The lowest breather is precisely degenerate with the soli- 
ton and antisoliton excitations, and they form a bound 
state which corresponds to S = 1 triplet. In this case, the 
elementary excitation should be a spin triplet and a spin 
singlet, no spinous exist as elementary excitation.^^"^^ 
However, the alternating NNN interaction, (dx^) cos 2$, 
does not lift the degenerate phases <l> = ±7r/4, due to the 
existence of the pre- factor dx^- 

We now perform a perturbative renormalization group 
(RG) analysis to our model by following standard 
procedures. Introducing the reduced variables ys = 



(10) 
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■ , we obtain the RG equations 

-ylK'+ylK' 

(2 - AK)y3 + 

(1 - K)yi - AK%m 



(7) 
(8) 
(9) 



under a change of length scale a ae'^''. Here we define 
dl = in°d^. The RG equation for the spin velocity u 
is a consequence of the anisotropy of the gi-interaction 
in the classical 2D XY-model, i.e. its non-retarded but 
non-local character in the quantum field theory (3). 

For 2/1 = equations (9) and (10) would not ap- 
pear, thus u is not renormalized. The RG equations (7) 
and (8) with j/i = describe the symmetric spin lad- 
der system with a Kosterlitz-Thouless transition.^ Spin- 
rotation invariant models scale along the separatrix be- 
tween Ncel and spin liquid phases, or along its con- 
tinuation into the dimer regime. Linearizing the RG 
equations (7)-(10) around the isotropic Heisenberg fixed 
point {g^"^ = 0, g["^ = 0, X(^) = 1/2) and defining 



6K = K 



2' 



we obtain the linearized RG equations 



dSK 

dl 

d ys 
dl 

d yi 
dl 
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-4SKy3 + -yi^ , 



(11) 
(12) 

(13) 



A family of solutions of the RG equations, projected 
on the 2/3 — /C-planc, are shown in Fig. 2. We choose 
the initial value of j/i as 0.001 and find the trajecto- 
ries are not sensitive to the choice of the initial values 
of 2/1. From the RG equations, one can directly find 
that there exist two intermediate fixed points given by 
{6K,ys,yi) = (0.1, 0.4, ±0.8). Here yi takes the values 
±, which reflects that our RG equations is symmetric to 
2/1. As shown in Fig. 2, the intermediate fixed point 
{5K*,y^) = (0.1,0.4) on the plane of 2/3 - SK is stable 
along the line 2/3 = 4:SK where the spin-rotation invari- 
ance is protected. For points near the intermediate fixed 
points, but not exactly on the line of 2/3 = 4:SK, the 
spin-rotation invariance is broken, thus they will flow to 
the spin-fluid fixed point (2/3 = 0) or the strong coupling 
fixed point (2/3 00). 
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FIG. 2. The scaling trajectories for yi{l = 0) = 0.001 pro- 
jected on the j/3 — 6K plane. 5K = K — 1/2, and the dot 
locates the new intermediate coupling fixed point. The Ncel 
state is realized in the upper left, the dimer state in the lower 
left, and the spin liquid in the right part of the figure. The 
asterisk locates the boundary between flows to the new fixed 
point, and into the dimer regime. 

At the fixed point, the RG equation (10) implies that 
the spin velocity u* is renormahzed. The robustness 
of the existence of this intermediate fixed point against 
higher order perturbations does not depend on the ex- 
act value of the fixed point as long as the fixed point 
is located on the RG scparatrix with 1/2 < K < oo.^ 
The intermediate fixed point is thus described by an ef- 
fective fixed-point Hamiltonian. Inverting the definitions 
of Ui = gi/iTu, we rewrite the effective fixed-point Hamil- 
tonian as a product of u* and a term independent of u*, 
H* = u*H{K*,g^,gi). Then the vanishing of u* leads 
to a trivial fixed point Hamiltonian H* = 0. As a re- 
sult of the vanishing of the renormalized spin velocity, 
the elementary excitations at the fixed point, spinon and 
antispinon, are still gapless. The vanishing of velocity is 
usually interpreted as a sign of ferromagnetism, however 
we interpret this as our spins effectively decoupling at the 
lowest energy scales, i.e. a kind of asymptotic freedom 
in this spin-rotation invariant ladder^. Also, the numer- 
ical results of Wiessner et al.^^ indicate a paramagnetic 
susceptibility. 

When J2 increases beyond a critical value J2c{^) (now 
depending on S) , the RG fiows to a strong coupling fixed 
point, which corresponds to the quantum dimer phase. 
For ui = 0.001, this critical point is indicated in {K,y^)- 
coordinates in Fig. 2 by an asterisk. For small 5 and 
J2 > J2c{S), our RG equations show that the system will 
remain in the universality class of the dimer solid cor- 
responding to the strongly fixed point, however the spin 
gap is decreased by increasing S. Basically, the gap size 
A cx exp(— Zi) where li is the length of the scaling trajec- 
tory from the initial values to the point where the most 
relevant perturbation is of order unity. This length is 
increased, and A therefore decreased, by the yi contri- 
butions to K and ys being opposite in sign to those of K 
and ys in Eqs. (7) and (8). 

Recently, Sarkar and Scn^^ studied the same model 
by using a nonlinear cr-model field theory and Abelian 
bosonization. However, the main discrepancy between 
our work and theirs^^ is that we kept the bosonized oper- 
ator of the alternating NNN operator and analyzed it by 
RG, while they just discarded it by giving an argument 
of the irrelevance of the operator. For an anisotropic 
XXZ chain, our RG result indeed shows that the gi op- 
erator is irrelevant^ in the meaning that it does not drive 
the system to a new phase, and this is consistent with 
that of Sarkar et. al. But the main difference lies in the 
question whether an intermediate fixed point exists and 
whether this fixed point corresponds to a phase differ- 
ent from a Luttinger liquid. In our previous work the 



argument of a vanishing spin-wave velocity plays a cru- 
cial role in the existence of such an unusual phase. If 
the spin-velocity does not vanish, one should explain the 
fixed point as a spin-liquid phase as in the Ji — J2 model 
with the spin velocity renormalized. Our result also 
suggests that the quantum phase transition parameter 
{J2c/Ji — 0.2412) to the dimer phase is changed by the 
alternating NNN operator, which may be verified directly 
by numerical simulations like the density-matrix renor- 
malization group. We also notice that the magnetization 
curve of the Heisenberg model with an additional alter- 
nating NNN operator"'^* gives an obvious different mag- 
netization susceptibility from the one without it. Their 
susceptibility is enhanced by this new interaction with 
respect to an equivalent Heisenberg chain, indicating a 
reduced spin velocity. Furthermore, a complete scheme 
to deal with the alternating NNN operator should give a 
correct description of the induced effect not only on the 
weak frustration regime but also on the strong frustra- 
tion regime. The omission of the gi term could not give 
any explanation why the operator shrinks the spin gap 
sizes in the regime of strong frustration as we will study 
in detail in the next section. However, our RG analy- 
sis gives a qualitative explanation of the infiuencc of the 
gi term on spin gap sizes. Therefore, we think that the 
scheme of discarding the gi operator based on its irrel- 
evance seems to be oversimplified. We hope that more 
numerical simulations will eventually be able to resolve 
this disagreement and study quantitatively the phase di- 
agram of the ground state as a function of J2/J1 and 
5. 



III. CROSSOVER FROM M-G TO SAWTOOTH 
MODEL 

It is generally believed that continuum field theory 
cannot give a good description for the behavior of the 
system far away from the critical point. In the case of 
J2 = O.5J1, the correlations extend only to a distance of 
one lattice spacing, thus the continuum field description 
is not a good approach. As we have shown, the phase 
corresponding to large J2 is the dimer phase. This is 
consistent with our knowledge from the models with ex- 
actly known ground states, say, Majumdar-Ghosh (M-G) 
model and sawtooth model. In this section, we prove that 
there exists a continuous manifold of Hamiltonians with 
dimer product ground states as long as J2 = Ji/2. 

We start with the asymmetric ladder model 



2JV 



1=1 



S/-S 



1+2 



(14) 



The M-G model and sawtooth model are corresponding 
to S = and 6 = J/2 respectively. 

For the M-G model^", the two linearly independent 
ground states, say, the left and right dimer ground states, 
are products of nearest-neighbor singlets, respectively 
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\^l)= UilJ + l], \^r)= + (15) 

I— odd I— even 

where [i,j] = {ai(3j — f3iaj)/\/2 denotes the singlet com- 
bination of spin i and j with the direction of dimers de- 
fined as i ^ j. Here q,; represents the up-spin and f3i 
the down-spin state at site i. \ ^l,r) also represent the 
degenerate ground states of the sawtooth model For 
the asymmetric ladder model, we notice that the NNN 
exchange alternation does not modify the product states 
of nearest-neighbor singlets 



Hs I ^l,r) 



^(-1)'(5 Si -81+2 



^l,r) = 0. (16) 



This is induced by the fact that the alternating NNN cou- 
plings along the upper leg and the lower leg of the ladder 
cancel out each other, when they operate on | $l,_r). It 
is obvious that | ^l) and | are eigenstates of the 
Hamiltonian (14). In fact, as we will prove, they are 
exactly the ground states of (14). 

To see this more clearly, we rewrite the asymmetric 
ladder model as a sum of projection operators Pi 



3/2. 



N 



3, J 



3, J 



2'2 



2'2 



,3/2 

2n 



with 



p3/2 ^ 



{Si-i+Si + Si+if 



■ (17) 



(18) 



Here, we introduce I to indicate the center position of 

three neighboring sites (^ — 1, I, I + 1). Such an opera- 
tor is a special case of the general positive semidefinite 
Lowdin's projection operators^'^ 



= n 



(Si 



Smf-S{S+l) 



s=s„ 



^maxi^max + 1) ~ S{S + 1) 



(19) 



where Smax and Smm are the maximum and minimum 
values of the total spin 5. 

As long as \S\ < ^, the coefficients in Eq. (17) are 
non-negative. Therefore, the Hamiltonian (14) is a lin- 
ear combination of projection operators with positive co- 
efficients. Since -Pf projects a state composed of three 
spins, {Si-i,Si, Si+i), into a subspace of total spin |, its 
eigenvalues are (if the total spin is i) and 1 (if the to- 
tal spin is |). By virtue of the properties of the positive 
semidefinite projection operator, whose lowest eigenvalue 
is zero, the ground state of the Hamiltonian (17) can be 
constructed by choosing states with such configurations 
that each projection operator has the lowest eigenvalues 

when operating on these states. It is easy to prove that 

1 ^l,r) are the exact ground states of the asymmetric 
ladder model (14). The ground state energy is indepen- 
dent of 6 and given by 



E„ 



(20) 



There is thus an entire manifold of Hamiltonians with 
fixed Ji = 2J2, parameterized by 6, with doubly degen- 
erate ground states of NN-dimer product ground states 
^L,R- For convenience, in our following discussion we will 

shift our model by an energy of Eg, H — Eg ^ H , which 
is equivalent in taking the ground state of the system as 
zero. 

In the following, we consider the excited state of our 
system (14). The elementary excitation of the system 
is a pair of spinous known as the kink or antikink''"^ . 
First results were obtained by Shastry and Sutherland 
for the M-G modeP^ within a variational ansatz. Since 
the ground state of (14) is independent of the alternating 
NNN exchange, the construction of the excited states for 
the M-G model can be directly extended to the asymmet- 
ric model (14). Breaking a singlet pair in the groimd state 
would give rise to two unpaired "defect" spins. There- 
fore, the simplest excitation consists of a pair of spinous. 
The spinons can be thought as domain-walls separating 
different dimer ground state configurations. From sym- 
metry consideration, the kink and antikink are identical 
in the M-G model, and specifically they have the same 
dispersions. With alternating NNN interaction, the sym- 
metry between legs is broken, therefore some properties 
of kinks and antikinks are different, in particular their 
dispersion. However, they still survive as elementary ex- 
citations of the asymmetric spin ladder system. 

antikink 
2n 

Xl^inRTS 





2m-1 



2n-1 



FIG. 3. The kink and antikink excitations in the asymmet- 
ric ladder model. The double lines represent singlets. 

In general, we call a spinon at the odd site 2m — 1 a 
kink and the other one at the even site 2n an antikink (see 
Fig. 3). The kink and antikink always appear in pairs 
in a periodic system, however, a single spinon can be 
realized in the open boundary systems. Taking the vari- 
ational wavefunction with one "defect" spin,^^ one can 
easily obtain the spinon dispersion of the M-G model 



5 J 
e{k) = - J + —cos2k. 
8 2 



(21) 



The energy gap A is therefore J/8. For S = J/2, i.e. 
the sawtooth model, the kink excitation is much differ- 
ent from the antikink excitation. As shown by D. Sen et. 
al. and T. Nakamura et. al.^'^, the kink (K) excitation 
in the sawtooth chain is exactly a single spin on odd site 
and dispersionless 



£K{k) = 0. 



(22) 
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However, an antikink propagates with an effective mass 
along the lattice. The antikink is not a free spin and 
spreads out to an extended region, because it is not a 
eigenstate of the local Haniiltonian. In the first approx- 
imation, the antikink (K) is supposed to be a single de- 
fect spin at the even site, and the dispersion obtained 
by variational calculation has similar form as the spinon 
dispersion of the M-G model 



£j({k) = -^J + Jcos2k. 



(23) 



with the corresponding energy gap of J/4. Despite the 
variational nature of the dispersion, the results agree very 
well with exact numerical results'^^^'^^. 

We will explicitly calculate the change of the gap size 
with increasing S. As expected, we found that the S 
term changes the energy gap size of excitations, which 
is consistent with our conclusion obtained by the renor- 
malization group analysis. Following D. Sen et. aF, we 
assume both the kink and antikink to be a 5-cluster block 
with spin 1/2. It is known that for the M-G and saw- 
tooth chain there is no closely bound kink-antikink pair 
whose energy is lower than that of a widely separated 
pair. Thus we can deal with the kink and antikink sepa- 
rately. The gap of the lowest excitation is a sum of the 
gaps of kink and antikink 



A = Ak + A- 



'■K ' 



(24) 



where the subscripts, K and K, represent the kink and 

the antikink respectively. It should be noticed that both 
the M-G and sawtooth model have the same energy gap 
size A under the first approximation (1-cluster approx- 
imation). That is not true as we take more precise n- 
cluster approximation. 

Under the 5-cluster approximation, the only three lin- 
early independent configurations that we need to consider 
are those shown in Fig. 4. We denote these three config- 
urations of kink by | 2to— 1 )i , | 2m— 1 )2 and | 2to— 1 )3 
and the configurations of antikink by | 2n )i , | 2n )2 and 
I 2n )3 respectively. Here, (2m — 1) and 2n denote the 
position of the center of 5-spin cluster corresponding to 
the kink and antikink. We now consider the momentum 
wavefunction with two variational parameters ai,2 and 

|fci) = -^^ gi(2m-i)fci [|2m-l)i+ai|2m-l)2 

^ m 

+ &i|2m-l)3] (25) 
1^2) = ^ E [|2n)r+a2|2n)2 



+ h2\ 2n )3 ] , 



(26) 



where ki and ^2 are the momentum of the kink and an- 
tikink respectively. 



2n 
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FIG. 4. The five-size cluster of the kink and the antikink. 

The lowest energy is obtained by finding parameters 
which minimize the energy expectation 



e(fci,2) 



(fcl,2|g|fcl,2) 

(fci,2|A;i,2) 



(27) 



Since |2m — 1)2 and |2m — 1)3 ( |2n)2 and |2n)3 ) are 
symmetric about the site 2m — 1 (2n), there is no rea- 
son to discriminate between these configurations and it 
is reasonable to choose ai,2 = ^1,2 ■ 

The computation of Eq. (27) is straightforward al- 
though a little bit lengthy, we will not give the detail 
here but refer to the literature^. It is found that the 
minimum value of e(fci) occurs at fci = 7r/2 and is given 

by 



e{ki = 7r/2;ai) 



1 J/2 (1 + 4ai) - 5 
4 l-ai+ al/2 



and the minimum value of £(^2) is 



e{k2 = 7r/2;a2) 



1 J/2 (1 + 4al) + S 
4 l-a2 + ay 2 



(28) 



(29) 



where we take J = 1 for convenience. 

For any given value of 6, the excited gap Ax i^-x) of 
a kink (an antikink) corresponds to the minimum of Eq. 
(28) (Eq. (29)). Our result is shown in Fig. 5, which in- 
dicates that the elementary excitation gap decreases from 
0.234 in the M-G to 0.219 in the sawtooth with the in- 
crease of the coupling constant while the ground state 
energy is constant. In particular, for 6 = 1/2, i.e. the 
sawtooth chain, Eq. (28) reduces to 



e{ki = 7r/2;ai) = 



2al 



2 /o ' 



4 l-ai+ ai/2 



(30) 



which has the minimum Ak = for ai = 0, while Eq. 
(29) becomes 



e(fc2 = 7r/2;a2) 



1 l + 2a^ 



4 1 - 02 -h ay 2 



(31) 



which has a minimum A-j^ = 0.2192 at 02 = —0.2808. It 
is clear that the kink excitation is exactly dispersionless, 
while an antikink is still a domain wall propagating with 



6 



an effective mass'''**. For the M-G model, 5 = 0, Eq. (28) 
and Eq. (29) have the same form {Ak = Aj^), thus 



£(7r/2;a) 



1 1 + 4a^ 



. + a2/2 ■ 



(32) 



whose minimum value is 0.2344 . In these limits, the re- 
sults are consistent with the known results of the M-G 
and the sawtooth model, as well as our qualitative con- 
clusion obtained from the field theory. 




0.1 0.2 0.3 0.4 0.5 

FIG. 5. The elementary excitation energy A versus S at 
the M-G point J2 = 0.5 with A = Ak + A-^. 

If we introduce additionally an alternating NN ex- 
change to our asymmetric model (14), the degeneracy 
of I $z,) or I will be lifted and the singlets would be 
pinned along the stronger NN external dimer potential. 
The elementary excitations are no longer separated kinks 
and antikinks. In the presence of an external dimer po- 
tential, a kink and an antikink separated by a distance 
of / give rise to a region in the incorrect "ground state" , 
which effectively produces a confining potential between 
the kink and antikink. This potential is proportional 
to the distance of I, thus the kink and antikink can not 
escape from each other and behave analogously as quark- 
antiquark pair.^^^^** The kink-antikink boimd state cor- 
responds to a magnon with spin 1 . The interesting topic 
of how the confined spinous develop to magnons has been 
investigated by Uhrig et. al.^^ 



exchange interactions and on top and bottom lay- 
ers may have different strengths. The Hamiltonian of our 
model is given by 

JV,M 2 

i,j=l a=l 
N,M 

N,M 



where the superscripts a = 1,2 denote labels of the top 
and bottom layers. J± is the perpendicular inter-layer 
exchange interaction and Jd is the diagonal inter-layer 
exchange interaction. Here all the exchanges are taken 
to be positive. A similar model has been investigated in 
Ref.^*, where the layer model is a direct generalization 
of the Bose-Gayrcn ladder modeP*'. It is clear that every 
slice of the double layer net is just a ladder whose Hamil- 
tonian has the same form of Eq. (14). Thus we find that 
the ground state of the layer model is given by a product 
of all perpendicular singlet pairs 



M,N 



D 



when the condition 

Ji_ = 2Jd = 2(ji + J2) 



(34) 



(35) 



is fulfilled. A rigorous proof can be made directly by 

representing the layer model as a sum of the projection 
operators as in the spin ladder case. The corresponding 
ground state energy is 



En 



--NMJ±. 
4 



(36) 



IV. DOUBLE LAYER MODEL 



Recently, it was found that the 2D Shastry-Sutherland 
model^^'^^ can be used to explain the experimentally re- 
alized material SrCu2{B03)2, hence such kind of model 
with exact dimer ground state^^ attracted much atten- 
tion again. In this section, we will show that the asym- 
metric spin ladder model can be generalized to a double 
layer model, whose ground state is a simple direct prod- 
uct of singlet dimers. 

The double layer model is constructed from two cou- 
pled spin layers shown in Fig. 6, where each layer has 
N X M sites and couples to the other layer by the inter- 
layer exchange interactions Jj_ and Jd. The intra-layer 




FIG. 6. The two-dimensional double layer model. 
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It is obvious that the dimer product state remains 
the groimd state for J± > 2Jrf. The properties of the 
ground state are independent of the specific values of 
and as long as the constraint condition (35) is satis- 
fied. Since the dimerized ground state is not degenerate, 
the lowest excitation is expected to be a triplet excita- 
tion, corresponding to breaking of a singlet bond, with a 
gap size proportional to J± . However, the many-particle 
excitation spectra might be very complicated because of 
the effective interactions among the triplet excitations^". 



V. CONCLUSIONS 

The spin-isotropic, asymmetric zig-zag ladders are 
studied using the field theory method and the varia- 
tional approach. When the leg exchange integrals are 
small compared with the NN exchange, the spin model 
is mapped to a revised double frequency sine-Gorden 
model. Renormalization group analysis shows that there 
are two fixed points, say, an intermediate-coupling fixed 
point and a strong coupling fixed point. In the weak frus- 
tration limit, the system is described by the intermediate- 
coupling fixed point with gapless excitations. The van- 
ishing of spin velocity at the intermediate-coupling fixed 
point is likely to indicate a decoupling of spins at low en- 
ergy scales. Apart from the isotropic separatrix, we find 
gapless spin liquid and gapped Neel states with easy- 
plane and easy-axis anisotropy. For large frustration, a 
more usual dimer solid phase is realized corresponding 
to the strong coupling fixed point. The RG analysis also 
predicts that the spin gap is decreased by increasing the 
leg-asymmetry 6. A continuous manifold of Hamiltoni- 
ans with the same singlet product ground state interpo- 
lates between the Majumdar-Ghosh model and the saw- 
tooth spin chain. Starting from the exact ground state 
wavefunction, we construct the variational wavefunction 
of the excited state and investigate the change of spin 
gap with the change of leg-asymmetry S. In the spirit of 
constructing Hamiltonian in the form of a sum of pos- 
itive semidefinite projection operators, extension to the 
double layer model is carried out. We propose an exactly 
solved two-dimensional double layer model with a ground 
state of a product of inter-layer dimers. 
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